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Abstract 

Universal Extra Dimension (UED) models tend to favor a distinctively heavier Higgs 
mass than in the Standard Model (SM) and its supersymmetric extensions when the 
Kaluza-Klein (KK) scale is not much higher than the electroweak one, which we call 
the weak scale UED, in order to cancel the KK top contributions to the T-paramcter. 
Such a heavy Higgs, whose production through the gluon fusion process is enhanced by 
the KK top loops, is fairly model independent prediction of the weak scale UED models 
regardless of the brane-localizcd mass structure at the ultraviolet cutoff scale. We study its 
cleanest possible signature, the Higgs decay into a Z boson pair and subsequently into four 
electrons and/or muons, in which all the four- momenta of the final states can be measured 
and both the Z boson masses can be checked. We have studied the Higgs mass 500 GeV 
(and also 700 GeV with ^/s = 14TeV) and have found that we can observe significant 
resonance with the integrated luminosity 10fb~^ for six dimensional UED models. 
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1 Introduction 



The Universal Extra Dimension (UED) scenario [Ij, in which all the Standard Model (SM) 
fields propagate in bulk of compactified extra dimension(s), is an attractive possibility whose 
simplest five dimensional realization on orbifold 5*^/^2, the minimal UED model (mUED), 
may account for the existence of the dark matter as the Lightest Kaluza-Klein Particle 
(LKP) |2j and can give a loose gauge coupling unification at around 30TeV See also 
Refs. [H E| for review on mUED. For the mUED, the latest analysis including the effects 
from second KK resonances gives the preferred Kaluza-Klein (KK) scale at around Mkk ~ 
1.3 TeV [S]. As is mentioned in this result strongly depends on the brane-localized mass 

structure, which is assumed to be vanishing at the UV cutoff scale [7j in mUED. 

One of the most important signature to establish the model would be the direct search 
of KK resonances at the CERN Large Hadron Colhder (LHC). See Refs. I8]-[I1] for mUED 
and Refs. |15]-[21| for 6D UED models. We note that some of them also pertain to the 
International Linear Collider (ILC), see also Refs. j22j-|27j. The LHC already puts a lower 
bound on the KK scale for mUED as Mkk > 500 GeV at the 95% CL from Mt2 analysis 
of cascade decay of first KK particles into the LKP |14^ I28j . It is noted that bounds on 
mUED [29] and T'^/Z^ UED [SD] from 6 57 processes claim Mkk > 600 GeV and 650 GeV, 
respectively. Again all the above bounds strongly depend on the KK mass splitting and flavor 
mixing patterns and hence on the boundary mass structure]^ In particular, we cannot see 
a decay product unless there are enough mass splitting among the first KK modes so that 
it becomes sufficiently energetic. In this paper, we present a complementary signal that is 
insensitive to such detailed boundary structure. 

In the SM, the electroweak data constrain the Higgs mass to be Mh ^ 170 GeV at the 
95% CL, see e.g. |321 133j . On the contrary, mUED prefer heavier Higgs when the KK scale 
is not much higher than the electroweak scale f ew — 246 GeV, namely, the KK scale should 
be Mkk > 800 GeV (300 GeV < Mkk < 400 GeV) at the 95% CL for Mh = 115 (700) GeV 
[3ll[35ll36] . We note that this is fairly model independent feature of a general UED model since 
KK top modes always contribute positively to the T-parameter and such an effect requires a 
heavy Higgs in order to cancel these KK top contributions by the ordinary negative log Mh 
dependence. In this paper, we call such a natural UED model without big mass splitting 
among electroweak, Higgs and KK scales: Mh ~ Mkk ~ O(IO^) GeV, the weak scale UED 
model. Concretely, we will pick up the cases: Mh = 330, 500 and 700 GeV. To summarize, 
existence of a heavy Higgs is a model independent prediction of the weak scale UED models 
in contrast to the relic abundance of the dark matter and the cascade decay signature of the 
KK particles that are dependent on the detailed boundary mass structure at the Ultra- Violet 
(UV) cutoff scale of the higher-dimensional gauge theory. 

In five and six dimensional UED models, the Higgs production cross section via the gluon 
fusion process is enhanced by the KK top loops |371 |38l [39]. In this paper, we analyze its 
cleanest possible signature, the Higgs decay into a Z boson pair and subsequently into four 
electrons and/or muons, in which all the four momenta of the final states can be measured 
and both the Z boson masses can be checked. 

In Section [2j we review the relevant part of all the known 5D and 6D UED models to the 

^ A 95% CL bound on the KK scale Mkk > 961 GeV is put on a "UED" model, assuming existence of large 
additional extra dimensions compactified with radius of order eV~^, in which SM fields cannot propagate, so 
that the LKP decays into KK-gravitons 31 . In this paper, we do not assume such additional large extra 
dimensions. 
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Higgs production process via the gluon fusion through the KK top loops. Considered models 
are the 5D mUED model on /Z2 [1], the Dirichlet Higgs (DH) model on an interval |401l41j. 
the 6D T^-based models on T^/Zs [H, T^/{Z2 x Z'^) I42j, T'^/Za 031 M^P^ [35], and the 
6D 52-based models on Projective Sphere (PS) [36], 5^ (see Section |2!3.2| ), 57^2 @7]0 In 



Section |3j we present the concrete computation of the process. The cross sections for the 
DH, rV2'2, TV(^2 X Z'2), RP"^, and S"^ UED models are newly obtained. We also review 
the estimation of the UV cutoff scale in the T^-based geometry |48j and extend it to the 5^- 
based one. In Section [4j we show our numerical results. The last section is for summary and 
discussions. In Appendix [Aj we present the relevant Feynman rules for our computation in 
the DH model. In Appendix [B] we explain our estimation of the UV cutoff scale for 6D UED 
models, based on the Renormalization Group Equation (RGE) analysis in the renormalizable 
KK picture. In AppendixjCj we review the way to take into account the width in the amplitude 
and justify our approximation. 

2 Review on known 5D and 6D UED models 

In this section, we give a brief review on various UED models. Readers who are not interested 
in the details of these models may skip this section. In the first part of this section, we 
briefly review the 5D minimal UED model on S^/Z2 [1] and Dirichlet Higgs model on an 
interval [401 3T] • The remaining of the section is devoted to an overview of various types of 
6D UED models. 

2.1 5D UED models 

2.1.1 Minimal UED model on S^/Z2 

First we review the 5D UED model [1]. The matter contents of the model are the same as 
those of the SM, but they are living in the bulk of flat five-dimensional space, compactified 
on the orbifold /Z2. The action S is written down as 

S= f (fx dy [£buik + 5{y)i^o + 5{y - ttR)C^r] . (1) 

J J~ttR 

Usually when one say mUED model, it is implied that all the boundary masses are zero at 
the UV cutoff scale and are generated through radiative corrections j7j. Hereafter, when we 
call mUED model, we do not assume any boundary mass structure and concentrate on the 
signal that is independent of it. In particular, we do not include the constraints from the 
direct KK search |14l [2H] and from the relic abundance of the LKP [S] that are dependent on 
the KK mass splitting pattern. 

The Z2 twist conditions on the bulk SM fields are put as 





-y) 


= B^,{x,y), 


B^{x, -y) = 


-B5{x,y), 




-y) 


= W^,ix,y), 


m{x,-y) = 


-m{x,y) 




-y) 


= Qii{x,y), 


G5{x, -y) = 


-G5(.x,y), 



(2) 

^ In |46) the terminology "real projective plane" is employed for the compactified space, the sphere with 
its antipodal points being identified. In order to distinguish 46 from 45 , we call the former the Projective 
Sphere. 
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L{x,-y) = 7^L(x,y), 
Q{x,-y) = j^Q{x,y), 



E{x,-y) = -j^E{x,y), 
U{x, -y) = --/^U{x,y), 
D{x,-y) = -j^D{x,y), 



and 



^{x,-y) = $(a;,y), 



(3) 



(4) 



where x and y {= x^) denote four and extra dimensional coordinates, respectively We can 
see that the wanted zero modes remain after the twist Q-Q. There are fixed points of 
the Z2 orbifolding at y = 0, ttR. If the boundary Lagrangians at y = 0, ttR are equal at 
the UV cutoff scale, there remains an additional accidental symmetry under the reflection 
^ — y — )• ^ + y, called the KK parity, which ensures the stability of the LKP and makes it 
a dark matter candidate. 

The gauge and Yukawa interactions for the KK-top quarks, which we need for later cal- 
culation, are: 



Ckk top = -ig4s J2 ^] 



n=l 



(n) 



mt 

VEW 



(n) 



n=l 



sin 2a(") 
7^ cos 



—7^ COS 2a(") 
5 — 2a(") sin2a(") 



in) 



(5) 



where g^s = (7s/\/2vri? is a dimensionless 4D SU{?>)c coupling constant and iiew — 246 GeV 
is the 4D Higgs vacuum expectation value which appear after the KK expansion; 
shows zero- mode gluon (zero- mode physical Higgs) ; t'^^^ and 4"'' are mass eigenstates of n-th 
KK top quarks and each mixing angle a^^^ is determined to be cos 2a^^^ = "^'(n) / \J'^t + ™\n) ' 



sin2a('^) 
n-th KK top mass is 



ml + m^ 



with m 



in) 



n/R. Each KK state is twofold degenerate and 



mf + 



(n)- 



(6) 



KK tops give dominant contribution to the gluon fusion process due to their large Yukawa 
coupling to the Higgs. We note that 7^ is put in Eq. ([s]) merely to arrange the sign of both 
the KK masses positive. 



2.1.2 Dirichlet Higgs (DH) model 

Dirichlet Higgs model is defined on an interval: < y < ttR. The action S is as follows: 

S = jd^x j^"^ dy [£buik + <5(y)£o + Siy - 7tR)C^r] , (7) 

where R is a radius of the extra spacial direction. The structure of the bulk Lagrangian, 
covariant derivatives and field strength of gauge bosons are the same as that of the mUED 
model. There is no difference between the matter contents of this model and those of the 

^ We follow the metric and spinor conventions of [49j . 
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mUED model. As in the mUED model, we neglect the possible boundary interactions in this 
paper. The zero-mode sector of the UED on an interval becomes the same as that of the 
mUED on the orbifold jZi when we choose the boundary conditions for the SM degrees of 
freedom = W^, L^, Ql\ Er, Ur, Dr to be Neumann (at y = and ■kR): 

55^^(x, 0) = d^^^{x, ttR) = (8) 

and for other non- SM modes ^'^ = B5, W5; Lr, Qr; El, Ul, Dl to be Dirichlet: 

^'^(x,0) = ^'^(x,7ri?) = 0. (9) 

We note that mode functions with Dirichlet and Neumann boundary conditions are not or- 
thogonal to each other, unlike the orbifolding on S'^/Z2j^ Kinetic terms turn out to be 
diagonal even though the expansion is not orthonormal. We can explicitly check that the 
non-orthogonality does not lead to extra mixing for spinors even after the EWSB because 
non-orthogonal terms drop out due to the 4D-chirality. 

If we had put the Neumann condition on the Higgs we would get exactly the same zero- 
mode sector as in the mUED model on /Z2- In the Drichlet Higgs model on interval, the 
EWSB is caused by a non-zero Dirichlet boundary condition on the 5C/(2)vi/-doublet Higgs 
field m]. We assume that the KK-parity is respected by the boundary conditions on 
the Higgs field too. The advantage of the Dirichlet EWSB is that we do not need to assume 
the negative mass-squared in the bulk Lagrangian nor the quartic coupling which is a higher 
dimensional operator in 5D. Throughout this paper, we consider the minimal case: V($) = 0. 
We list the necessary Feynman rules in Appendix [A] 



2.2 6D UED models based on 

We consider a gauge theory on six- dimensional spacetime x T^, which is a direct product 
of the four- dimensional Minkowski spacetime and two-torus T^: < y < 27ri?y,0 < z < 
27rRz. We assume that the two radii of have the same value R = Ry = Rz for simplicityj^ 
When we use 6D Weyl spinor for 6D UED model construction, there is a constraint on the 
choice of 6D chiralities. The origin of this constraint is the cancellation of 6D gravitational and 
SU(2)i global anomalies that cannot be removed by use of the Green-Schwarz mechanism. 
This constraint requires the number of matter generation to be (multiple of) three |51j . A 
suitable choice of the 6D chirality for a single matter generation is as follows: 

{Q+,U-,D^;L+,E^,N^), (10) 

where the it suffixes represent 6D chirality of each field. Number of d.o.f. of 6D Weyl fermion 
is 4, the same as that of a 4D Dirac fermion. Therefore we can construct 6D UED models on 
following the orbifolding method of the 5D UED model. We have several options for the 



orbifolding to realize the SM chiral fermions in the zero mode sector of (10). Let us review 
them in turn. The range for KK summation is listed in Table [2j 
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type of orbifolding 


identification 


fixed points {yi,Zi) 




{y,z) ~ {-y,-z) 

{y, z) ~ {-y, z) and (y, z) ~ (y, -z) 

(y,^) ~ {-z,y) 


(0,0), (7ri?,0), (0,7r/?), {-kR.ttR) 
(0,0), (7ri?,0), (0,7ri?), {-kR,ttR) 
(0,0), {ttR^ttR) 



Table 1: Fixed points whicli stem from each identification. 



2.2.1 Orbifold rV^2, T'^/{Z2 x Z^), and T'^/Z^ 

, and TVZ4 



We consider rVZs [Ij, rV(Z2 x ^ 
the action 



orbifolds. Let us write down 



S 



ttR 



dy 



TVR 



ttR 



dz 



ttR 



>Cbulk(a:,y,2:) + 



e »7i 



5{yi)Cg^{x) 



(11) 



where yi = {yi, Zi) are orbifold fixed points. We note that terms localized at the fixed points 
are induced at quantum level even if we assume that they are vanishing at tree level [3 E21 E3] • 
For the orbifold we consider in this paper, the projections are: 

(y,z)~(-y,-z) forrV^a, (12) 

{y,z) ^ {-y,z) and {y,z) {y,-z) for r^/(Z2 x Z2), (13) 

{y,z)r.{-z,y) forrVZ4. (14) 

See also Table. [T| In each case, we can choose a suitable boundary condition of 6D Weyl 
fermions, whose exact forms are not discussed in this paper, to generate 4D Weyl fermions at 
the zero modes. 

The bulk Lagrangian, covariant derivatives and field strengths of gauge bosons are essen- 
tially the same as that of the 5D mUED model except for the structure of spinors. For 6D 
Weyl fermions, the kinetic and Yukawa terms are 

/:kinetic = -q:^t^'VmQ+ - TjZt^'VmU- - dZt^'VmD. 

- L^T^^VmL+ - EIT^^VmE- - NIT^'VmN^ , (15) 
^Yukawa = -^uTT (Q+ • $) - Afl (Q7^>) - Xe (L^^) E. + h.c, (16) 

where contraction of SU (2) indices are understoodj^ Resultant interactions relevant for our 
discussion are 

C-KK top ■ 



Wis ^ [tl t2 

(m,n) 



{m,n) 



{m,n) 



mt 

VEW 



{m,n) 



{in,n) 



sin2a('"'") 
7^ cos2a(™'") 



-7^ cos 



sin2a('^''^) 









*2. 



(m,n) 



(17) 



■* In other words the KK mass-squared operator 9| is not hermitian in this setup, though the kinetic term 
is stiU positive definite. 

^ In the T^/Zi orbifold case, the condition Ry = is imposed by the consistency with the Zi discrete 
symmetry. See also [SD] for a realization of CP violation from the complex structure of T'^ /Z^, which appears 
in 4D effective interactions after KK decomposition. 

^ We leave the neutrino sector untouched since it is irrelevant for the Higgs signal considered in this paper. 
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type of orbifolding 


range of (m, n) 




m + n > 1, or m = —n > 1 
< m < oo, < n < oo; (m, n) ^ (0, 0) 
1 < 771 < oo, < n < oo 



Table 2: The range of the parameter (m, n) except the zero mode case (m, n) = (0, 0) in each 
case of the orbifolding. 



where = gs/ (27ri?) is the dimensionless 4D SU{^)c coupling constant and uew — 246 GeV 
is the 4D Higgs vev, 

g(0)(-^(0)^ shows zero-mode gluon (zero-mode physical Higgs), and 
^(m,n) ^ ^{rn,n) xnass cigenstatcs of (m, n)-th KK top quarks. Again we only consider the 
KK top quark loops since contributions from other flavors are suppressed by the small Yukawa 
coupling. Each mixing angle a^*"'") is determined to be cos2a('"'") = ?n'(m,ra) / \J'^'t + ™\m n) ' 

sm ml + Each KK state is twofold degenerate and (m, n)-th KK top 



mass IS 



(m.n) ' 

with 



mt,^m,n) = \ + mf., (18) 



-|- 

m(m,n) ■= ^ • (19) 

It should be mentioned that the difference from the mUED case appears only in the form of 
KK mass and the number of d.o.f. in each KK level when we consider the gluon fusion process. 
We adopt m(n) as the y(z)-directional KK index, whose parameter region is determined by 
the way of the orbifolding. This information has a great influence on the enhancement of the 
Higgs production through the gluon fusion. 

2.2.2 Real Projective Plane (iJP^) 

We can construct a UED model on a non-orientable geometry: Real Projective Plane {RP^) [35] 
RP^ is defined by two types of identifications: a vr-rotation r and a glide g: 

r : {y,z) {-y,-z), g : {y, z) {y + ttR, -z + nR). (20) 

The system is invariant under each manipulation in Eq. ( |20[ ) . Note that the shifts y ^ y + 2ttR 
and z ~ z + 2TrR can be obtained as different combinations of r and g, respectively. Note also 
that no fixed point exists globally in this background geometry. 
Under r and g, Weyl fermions transforms as 

r:^±{x;-y,-z)=prrr^±{x;y,z), = ir^r^r^ (21) 
g:^±{x;y + 7rR,-z + TTR)=pgTg^^{x-y,z), = T^r^ (22) 

where F'' is the 6D chirality operator and pr,Pg (Z2-parities) can take the value ±1. The ^± is 



what we call the "mirror" fermion. Eq. (21 ) has the same form with that of the T /Z2 orbifold 



condition for 6D fermion. An essential point of this model is that the condition ( 22 ) does not 
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generate a 4D Weyl fermion in the zero mode sector. In other words, the 6D chirahty of both 



sides of Eq. (22) are different from each other. This means that we have to introduce new 



fermions which have opposite 6D chirahty and the same SM quantum number compared 
to each corresponding field ^z^. Concretely, "mirror" fermions: 

Q.,U+,V+;L_,£+,M+, (23) 

are identified with {Q-^-,U-, D^; L-^-, E^, N^}, respectively. The choice of 6D chiralities in 



Eq. ( 23 ) obeys the condition for realizing the 6D anomaly cancellation which we have argued 



before. 

The bulk Lagrangian is the same as that of the T^-based models using orbifold except for 
the existence of the mirror fermions: 



inetic 



Yukawa 



- Ad (QT^) V+ - Xe (lI$) £+ + h.c. 



(24) 



(25) 



where we introduce the "1/2" factors for later convenience. The neutrino sector is again 



left untouched as it is irrelevant for our discussion. By use of Eq. (22), we can erase all the 



mirror fermions and obtain the ordinary form of Lagrangian same as cases. The form of 



vr-rotation r given in Eq. (21) is the same as that of Z2 orbifolding in Eq. (12). Therefore, 
the interactions of RP^ model, needed to calculate the gluon fusion process, take the same 



form as that of T /Z2 one given in Eq. (17). 



2.3 6D UED models based on 

Let us review UED models based on the 5^ compactification. We span the extra dimension 
by the zenith and azimuthal angles 9 and (p, respectively. The two-sphere 5^ has a positive 
curvature and to stabilize the radius R, we introduce an extra U{l)x gauge field which has 
a monopole-like classical configuration 



n 



2gx 



(cos 6* =F 1), 



(other components) = 0, 



(26) 



where the superscript c denotes the classical configuration, gx is the 6D U{l)x gauge coupling, 
the integer n is the (negative) monopole charge, and the superscripts N and S indicate that 
the field is given in north (involving the 9 = point) and south (involving the 9 = tt point) 
charts, respectively. The U{l)x transition function from the north to the south chart is given 

by 

[A'm(x^ 9, <p)f = [A'M(x^ 9, cP)f + — 5A/a(x^ 9, </.) (27) 

9x 
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with a{x^,9,(j)) = ncj). Because of the monopole-hke configuration, the radius of 5^ is stabi- 
hzed spontaneously at 



where M^, is the 6D Planck scale. 

We mention that any 6D field H on is KK expanded by use of the spin-weighted 
spherical harmonics sYj„i(^, (p) as follows: 

EixA<t>fs = f:f: f^'-\e,^fs ■.= sY,M)e^-\ (29) 

j=\s\ m=-j 

where is the (j, m)-th expanded 4D field, f!^'"^^ is the corresponding mode function and 

s is the spin weight of the field H. The spin-weighted spherical harmonics sYjm(0, (p) matches 
the orthonormal condition as 

J d(t) j dcOsO sYjm{0,(j)) sYfm'{0,4>) = djf5mni'- (30) 

A spin weight of fermion is closely related to its U{l)x charge. When we assign U{l)x 
charges of 6D Weyl fermions ^± as q\f/_^ , the corresponding spin weights of 4D Weyl fermions 
{■0 I L , ip_ L } are given as follows in our convention: 

.^.=^^, «_^=^^. (31) 

Note that if a 6D Weyl fermion takes a spin weight s = 0, a j = mode appears as a 4D Weyl 
fermion with vanishing KK mass. This means that we can get chiral SM fermions without 
orbifolding in the case of S*^. When we take the values: 

(s+i?, s+L, S-R, S-l) = (0, 1, 1, 0), (32) 

we can create the same situation as in the T^-based models discussed before. The spin weight 
of the 4D-vector component of a 6D gauge boson is always s = and then there is a zero 
mode that can be identified as the SM gauge boson. On the other hand, extra dimensional 
components of the 6D gauge boson are expanded by the \s\ = 1 spin- weighted spherical 
harmonics and has no zero-mode. 

In our configuration, any (j, m)-th KK mode has the KK mass: 

= Ji • (33) 

An important point is that the form of the above KK mass is independent of the index of m. 
This means that there are 2j + 1 degenerate modes for each j. Note that the lightest KK 
mode has the mass y/2/R. 

As discussed above, the 4D-vector component of a 6D gauge boson has a zero mode. This 
is the case for the extra U{\)x gauge boson too. Phenomenologically the existence of an extra 
f/(l) interaction, under which SM fields are charged, is problematic |46) . In the following, let 
us see how to get rid of this massless U{l)x vector. 
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2.3.1 Projective Sphere (PS) 

We can construct a UED model compactified on the Projective Sphere (PS), a sphere S"^ 
with its antipodal points being identified by {9, 0) ~ (vr — 0, + tt) [46J. In the UED model 
based on PS, the 6D action takes a different form from that of the 6D orbifold UED models. 
One of the remarkable points of this model is that there is no fixed point on the background 
geometry PS. As in the RP^ model, we introduce "mirror" 6D Weyl fermions: 

Q_,ZY+,P+;L_,£:+,AA+ (34) 

which have opposite 6D chirality and opposite SM and U (l)x charges, compared to the fields 
{Q-\- , U- , -D_ ; , E- , N- } . Because of the existence of mirror fermions the kinetic term takes 



the same form as in the RP^ model ( 24 ) and the Yukawa interaction is modified to 
1 



■^Yukawa — 2 



- XuU^ (Q+ ■'^)-Xd {Q+^) - \e (L+$) 

- XyU^ (Q_ ■^)-\*D (2^^) ^+ - (LT^) + h-c] . (35) 



Like the RP^ case which we have discussed before, we introduce the "1/2" factors for a later 
convenience. The covariant derivatives in this model are given as 

Vm = dM + igsQljT^ + igmiT'' + igvBMY 

(for $), (36) 
Vm = dM + igsQliT^ + igy^liT" + igyBMY + igxq^i^l^ + ^m) + 

(for Q+,C/_,L>_;L+,^_,iV_), (37) 

Vm = Bm + igsQMl-T^^ + i5W^[-r"]T + igyBMl-Y] + igxq^{Xt, + Xm) + ^m 

(for Q_,^+,P+;L_,f+,AA+), (38) 

where VIm is the spin connection. The covariant derivative of Higgs is the same as that in the 
S"^ /Z2 case, but there is a difference between those of fermions and these "mirror" fermions. 
We discuss these points shortly below. 

As we mentioned before, projective sphere is a non-orientable manifold and has no fixed 
point. Let us consider the 6D P and CP transformations. Under the antipodal projection, 

' A'^(x,7r-0,0 + 7r)5 = X(^{x,e,(l))N ^ 

N S 

Xe{x,7T-e,(l) + TT)s = -Xl^{x,e,(t>)N, (39) 

N S 

{XI,X^}{x,7r-e,(^ + 7r)s = {{X^f ,X^}{x,9,<f>)N , 

where the superscript C denotes the 6D C transformation. Recall that the superscript c 
denotes the classical configuration. These conditions leave the monopole-like configuration 
invariant under the antipodal identification and projects out the unwanted U{l)x 4D-vector 
zero mode. In contrast, identification of a SM gauge boson A^^j should be done by another 
condition since we want the corresponding 4D-vector zero mode, where i shows the type of 
gauge group. We adopt the 6D P transformation and those identifications are written as 

' A^\x,7r-e,^ + TTfs = Ali\x,9,cf>)N, 
Af{x,Tr-9,(t> + 7r)s = -Af {x,e,(t>)N , (40) 
Af {x,7r- 9,(1) + irfs = ^»(x,e,0)^, 
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where it is evident that the zero mode of A^ff^ survives. We also identify Higgs with the 6D 
P transformation to obtain its zero mode: 

JV s 

<^{x, TT - 9,(f) + Tr)s = 9,(p)N. (41) 

Finahy, we discuss the identification of 6D Weyl fermions. Since 6D Weyl fermions have 
U{l)x charge and interact with the U{l)x gauge boson, they should be identified by the 6D 
CP transformation. The specific form of the 6D CP transformation, for example in the case 
of U- , is as follows: 

U+{x,7: -6,()) + 7r)'^ = PUS{x,e,^)N^ (42) 
where the matter field U- is identified to the mirror We decide the forms of covari- 



ant derivatives (37) and (38) on the criterion of invariance of the action under the 6D CP 



transformation. Using the identification conditions (39)^(42), we can see that the mirror 
fermions drop out of the action after the identifications and eventually we obtain the usual 
type of UED model action. This can be interpreted that all the modes of mirror fermions 
{Q_,Z^+, P+; L_, A/+} are erased and no mode of Z)_; L+, A^_} is projected 

out. The interaction terms which we need for calculation in this model is the same as those 



in Eq. (17). Only difference is the number of degenerate top KK modes in each j-level. 



2.3.2 52 UED with a Stueckelberg Field (5^) 

As a solution to the massless U{l)x problem, we can simply give a Stueckelberg mass [55), [5^, 
see also [S71ISH] for reviews, to the U{l)x field. We can make the unwanted U{l)x 4D-vector 



zero mode to be massive while preserving the classical monopole structure in Eq. (26). This 
way, we can formulate a UED model on S"^ with no field identification. Let us call this simple 
model the S'^ UED model. In the S'^ UED model, the matter contents, bulk Lagrangian, 
definition of field strengths and covariant derivatives and the configuration of the classical 
U{l)x field are the same as the PS model after removing the mirror fermions, except for the 
Stueckelberg field part. In contrast to the S'^ jZ^ orbifold below, there are no fixed point nor 
a localized Lagrangian anywhere on 5"^, as in the case of PS UED. The d.o.f. of KK fermions 
has no difference between the S"^ UED and PS one (after the antipodal projection). There 
is no need for an additional computation; all we have to do is to borrow the PS result as a 
whole when we are only interested in the gluon fusion process. 

2.3.3 5V^2 orbifold 

Although the above S"^ UED model with a Stueckelberg U{\)x mass is already phenomeno- 
logically viable, we may further perform a Z2 orbifolding on it [47] On this orbifold, the 
point (0, i;^) is identified with (tt — ^, —</>). The 6D action S is as follows: 



/■27r 

S = \ d^x \ dQ 



JO 



(43) 



^ This extra Z2 cannot project out the U{l)x gauge field. In ^7], the U{l)x is assumed to be broken by 
an anomaly. Since we need a classical configuration of the f/(l)x, it would be theoretically preferable to break 
it by a tiny Stueckelberg mass. 
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where \/—g = -R^sin^. This system has two fixed points of the Z2 symmetry at (^,<?!)) = 
(|, 0), (f , tt) and we describe the locahzed terms with vC(7r/2,o)) ^{■k/2,-k)i respectively. Like the 
T^-case, we do not discuss those parts in this paper. 

We can easily construct mode functions of 5*^/^2 fg'^f"\9,(j)) with spin weight s in both 
north and south charts following the general prescription [59] as follows: 

^ [sYjmie, (I)) + {-iy-%Y,_m{e, </>)] e±^^<^ for t = +1 
T , (44) 

^ ^ [sYjm{e, cP) - {-iy-%Yj.m{e, cP)] e±^^<^ for t = -1 



where t = ±1 is the Z2 parity. These mode functions have the property that /^"'^'!!!].x(vr — 

9, —<p) s = ±/j"'j'™2i(^) ^ • To realize the Z2 symmetry, we identify a field at {9, (f>) in north 
chart with the same field at (vr — 9, —(p) in south chart. 

The range of the summation over m shrinks from [—j, j] to [0, j] after the Z2 identification. 
Under the transformation of {9, (p) — t- {9,(p + tt), mode functions behave as 

fi^U,{0,cP+^fs = {-irfi!f^l^,{9,cpfs, fi^l_,i9,cp+nfs = -(-l)-/BS,=-i(^, 

(45) 

After some fields redefinition, we can find that each KK field has a KK parity (—1)™, which 
is a remnant of the KK angular momentum conservation. 

We focus on the m = modes of each j level. When we see the concrete forms of mode 
functions in m = 0, which are 

/Sri|i(^,0)^ =^(l + (-in-olSo(0,</'), (46) 



2R 

fit:T?=-ii0,^f' = ^(1 + • iYAO,^)e^"^, (47) 

/ir_\tLi(^,0)" = 2^(1 + (-1)^') • -iY,o{9,cP)e^'t (48) 

where we find that m = modes appear only in the case of even j. Then degeneracy of KK 
masses is 

j + 1 for j : even, 
j for j : odd, 

since m runs from to j. Again, this mode counting is the only important point when 
computing the enhancement of the Higgs production via gluon fusion process. We do not 
discuss the form of interactions which we need for calculating the gluon fusion process because 
there is essentially no difference from the T'^ case (17). 



Higgs production and decay into four leptons in UED mod- 
els 



In the SM, the cross section for the Leading Order (LO) one-loop Higgs production via the 
gluon fusion process and its subsequent decay into a Z boson pair: gg ^ H ^ ZZ is given 
by [60]. The LO parton-level cross section shown in Fig. [l]is: 



a 



4s 



2567r3 



mz 



1 + 



2mlf 



8m| 



+ A2 



4m| 



\I{s)\\ (50) 
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Figure 1: Feynman diagram which describes the dominant contribution to the gluon fusion 
Higgs production process and the subsequent decay to 2Z. 



where m]y,mz,mt, and Mh are respectively the W, Z, top quark, and Higgs boson masses, 
ctis = g'ls/^'^ is the 4:D QCD gauge couphng, s is the center-of-mass-energy-squared of the 
scattering partons, we employ the normalization for the Higgs vev: v\^^ = \l\piGp ~ 
(246 GeV)^, and the loop functions are defined as 



dx 



/(A) = -2A + A(1-4A) / =ln 

/(A) = A /"'—In 
Jo X 



x{x — 1) 
A 



+ 1 



x{x — 1) 
A 



+ 1 



(51) 
(52) 



Explicit result of the integral is 



dx 



In 



x(x — 1) 



+ l-ie 



1 



n 2 



arcsm 



In 



4AJ 
1 + V1-4A 



1 - VI - 4A 



(for \>\: 

(for \<\] 



(53) 



We have also defined / for later use for the Dirichlet Higgs model. In Eq. (50), we have taken 
into account the total decay width of the Higgs in its propagator: 



A = MhTh. 



(54) 



In the current analysis, we take into account the Higgs decay into VF, Z and top quark pairs, 
which are dominant when we consider the heavy SM Higgs boson: Mh > 2m\\r. Explicit form 



is shown in Eq. (124) in appendix. Note that we take into account only the top quark loop in 



the SM cross section ( 50 ) , given by the diagram shown in Fig. [T| since the Yukawa coupling 
to others are negligible compared to the top one. We have also ignored the contributions from 
the sub- leading box diagrams [61]. See Appendix [C] for further discussion on how to take into 
account the width. 



3.1 Gluon fusion process in UED models 

In this paper, we consider the KK-top loop contributions to the gluon fusion process in several 
UED models, namely, 5D UED model on 57^2 (mUED) [1], Dirichlet Higgs (DH) [IQ], 6D 
UED model on rVZa [T], T'^/Z^ ^ E], TV(Z2xZ^) [12], Real Projective Plane {RP"^) 
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Figure 2: A schematic description of the dominant contribution to the gluon fusion Higgs 
production process and the subsequent decay to 2Z. The hned blob indicates the effective 
vertex. 

g(0) g(0) g{0) 




g{0) g(0) g{») 



Figure 3: The effective vertex which describes the Higgs production from the gluon fusion. 



52 [47J, Projective Sphere (PS) |46j and S'^ with a Stueckelberg Field (S'^). We have given 
a brief review on these models in the previous section. The contribution from KK-top loops to 
the gluon fusion process is analogous to that of the top-loop in SM and the difference resides 
only in the loop function. Relevant Feynman diagram is shown in Figs. [2]and|3| The effective 
vertex, which is represented by the lined blob in the diagram, includes the contributions to 
the gluon fusion from the zero mode top quark and the KK top quarks. Since the zero mode 
sector of UED model regenerates the SM configuration, the result of the former contribution 



is the same as that of the SM in Eq. ( 50 ) . The forms of the latter contribution will be shown 
soon later. We note that the letter "iJ" in Figs. [2] and [3] shows the zero mode physical Higgs 
boson in all the cases except for the DH model where H stands for the first KK Higgs. 

For each model, we get the following result, where Jmodei indicates the corresponding loop 
function: 



;.model < (mz\ l+{s- 2m|) s_ J^_'^K\^ 

(55) 
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m+ri> 1 
or m— — n>l 



m>l,?i>0 



2 / m? 



m>0,n>0, 
(m,n)7^(0,0) 



s 



2 / m? 



s 



where the KK top masses are given by 



mt{n) '■=\l'^t+^ = \Jml + ^^^KK' (63) 
/ TTiP' -\- up' I 

"^i(m,n) := Y ml + — ^2 — = \Jml + {m?+ n"^) j^, (64) 



Here the Mkk is the first KK mass, which is written as 

Mkk = ^ (66) 

for the compactifications based on jZ^-, interval, and (namely, the mUED, DH, T'^ jZ^-, 
RP^ , T"^ /Z/i, and T'^ /{Z2 x Z'2) models) and is written as 

Mkk = ^ (67) 

for the 5^-based ones (namely, the 5^/Z2, PS and S"^ models). The gluon fusion process 
for mUED in IZ2 is first shown in Ref. [37] and for S"^ /Z2 in Ref. [38]. Also it has been 
calculated for T'^/Z^ and PS in Ref. [39]. The results for DH, TV^2, RP"^ , TV(^2 x Z^) and 
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5^ are newly presented in this paper. The factor \ple\ in Eq. (57) is equal to 2^/2/11 ~ 0.9. 
The origin of this suppression factor is non-orthonormality of mode functions on an interval. 
In the case of S'^-based compactification, there are some degenerated states, the number of 
which is described with n(j) on S"^ /Z2 and {2j + 1) on PS or S"^ for each KK-index j. The 



specific form of n(j) for the orbifold S /Z2 in Eq. (61) is as follows: 



ji + 1 for j: even, 
I J for J : odd. 

Several comments are in order. 

• The origin of the factor 2 in front of each KK summation is the fact there are both left 
and right handed (namely, vector-like) KK modes for each chiral quark zero mode. 

• All the KK contributions are positive and hence always enhance the Higgs production 
rate via the gluon fusion process, except for the DH model in which the zero mode Higgs 
contribution is absent. 

• Each value of Yukawa couplings of KK quarks to the Higgs is the same as that of the 
coupling between the corresponding zero mode fermion and the Higgs. We only consider 
triangle loop diagrams of the SM top quark and its KK excited modes because their 
Yukawa coupling to the Higgs is dominant compared to that of other fermions. 

• In each KK summation infinite numbers of KK modes contribute to the process in 
principle. In 6D, these summations are divergent and a suitable scheme of regularization 



is required, jmax in Eqs. (61) and (62) shows an upper bound of the summation over 



the index j. Further discussion will be shown in the following subsection. 



K in Eq. (55) is the so-called K-factor, a phenomenological approximation in order 
to naively take into account higher order QCD corrections. One may take K ^ 2 for 
Tevatron and K ~ 1.5-1.6 for LHC, respectively [32j. In the limit where the KK- 
loop is viewed as a contribution to the effective Higgs-gluon-gluon coupling, the QCD 
corrections to Higgs production are very similar between the SM and the new physics 
contributions. The reason is that the Higgs-gluon-gluon coupling always has the same 
structure, and only its coefficient changes. This is discussed in detail e.g. in Ref. |62| 
in the context of SUSY (but it works the same way in UED as in SUSY) . Therefore we 
have included a K factor also for the new physics terms as in Eq. ([55|). 



When the compactification radius is too large, namely, when the first KK W is lighter 
than half the Higgs mass, the Higgs can decay into a pair of KK particles and its decay 
width Th becomes broader. In this situation it becomes harder to find the evidence of 
the Higgs boson and hence we restrict ourselves within the region where the Higgs mass 
is smaller than twice the first KK W mass so that such a decay mode does not open up. 

Although the one-loop gluon fusion process is the dominant production channel of the 
Higgs, its contribution to the Higgs total decay width is smaller at least by three orders 
of magnitude comparing to the decay into a W pair in the case of the SM with Mh > 
300 GeV [32j . Even after enhancement of 0(10) from KK-top contributions, decay into 
gluon pair is still negligible. 
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• In Eq. (55), the Higgs decay width is taken into account by the naive Breit-Wignar 
formula in the denominator. When the Higgs mass is large, say Mh = 700 GeV, the 
Higgs decay width is as large as 180 GeV. In some literature, the expression Mh^h 
in the the Breit-Wignar formula is replaced by sTh/Mh. In Appendix [C| we discuss 
reliability of our treatment. 

3.2 UV cutoff scale in six dimensions 

In 6D UED models, since the gluon fusion process is UV divergent, we must consider upper 
limit of the summations of KK number in such models]^ First let us briefly review how the 
Naive Dimensional Analysis (NDA) is applied to the higher dimensional theory. Following 
the concept of NDA, a loop expansion parameter e in D-dimensional SU(N) gauge theory at 
a scale /i is obtained as 

1 27r^/2 

where Ng is a group index, g^i is a dimensionful gauge coupling in D-dimensions and A is a 
UV cutoff scale. The index i is introduced to express the type of gauge interaction and the 
remaining part originates from D-dimensional momentum loop integral. The cutoff scale A 
is the scale where the perturbation breaks down e(A) ~ 1. 

Precisely speaking, the dimensionful higher dimensional gauge coupling does not "run." 



Let us explain what is meant by the running coupling in (69), basically following Ref. [48]. 
When we consider a 6D theory (D = 6) with two compact spacial dimensions, an effective 4D 
gauge coupling g^i emerges after KK decomposition: g^ = g^i / \A^; where V2 is the volume of 
two extra dimensions. Concretely, V2 = (27ri?)^ and AttF? for and S"^, respectively. In this 
paper, we employ a bottom-up approach for the running gauge coupling. At energies below 
the first KK scale, theory is purely four dimensional (after integrating out all the massive 
modes of order KK scale) and the gauge coupling runs logarithmically. Let us then increase 
the energy scale. Every time we cross a KK mass scale, there open up the corresponding KK 
modes to run in the loops in the gauge boson two-point function. In all the scales, the theory 
is renormalizable and the running of the gauge coupling is logarithmic. However, due to the 
increase of the number of particles in the loops, the running of the gauge coupling becomes 
effectively power-law at the energy scales much above the first KK scale. This way, we get 
the effective power-law running of the gauge coupling, within purely renormalizable approach. 
In this paper, we neglect possible threshold corrections at the UV cutoff scale, since we are 
interested in what is the highest possible A that can be consistent with the low energy theory. 

In the above stated strategy, we get the following running of the 4D effective gauge coupling 
strength an{^) 



a 



uSM k6D , 6D 

.Mfi) ~ a^Mmz) - ^ In ^ + 2C7 ^ In - C ^ 



AfKK 



^ ' 1 



(70) 



where C = 7r/2 and 1 for T and S , respectively. In Eq. ( 70 ) , we have approximated that all 
the masses are degenerate in each KK level. Depending on models, some fraction of the KK 
modes are projected out, but we assume that all of them contribute to the running, in order 
to give the most conservative upper bound on the UV cutoff scale. More detailed explanation 
is given in Appendix [Bj 

* In 5D UED model, we can execute this mode summation with no divergence. 
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T^-based 
max mill 


S'^-based 
max mill 


KK index 
UV cutoff 


m'^ < 28 m'^ + < 10 
AgD ~ 5Mkk AgD ~ 3Mkk 


j(j + 1) < 90 j(j + 1) < 30 
AgD ~ 7Mkk AgD ~ 4Mkk 



Table 3: Our choices of maximum and minimum upper bounds for KK indices and for the 
corresponding UV cutoff scale. 



While the described procedure gives a reasonable estimate of the UV cutoff scale, one has 
to be aware that this is not much more than an order-of-magnitude estimate. We will plot 
our results for maximum and minimum values of the UV cutoff scale that are theoretically 
reasonable. In Table [3j we list our choice of bounds at which the KK mode summation is 
truncated to regularize the process. 



3.3 Convolution of parton distribution 

Let us briefly review the standard prescription to estimate the event number of the Higgs 
production in pp — )• ZZ — )• Ai (four leptons) via the gluon fusion process as a function of 
the invariant mass of ZZ, given the parton level cross section, where M denotes two pairs 
consisting of either e~e^ or fi^ fi^ , since a tau pair is less visible at the LHC. That is, the final 
state is possibly e~e~^e~e~^, fi^ fj,^ fj,~ fi'^ , or e~e"'"/i~^"'". Because of the large mass difference 
between Z boson and two electrons or two muons, the subsequent processes: Z e+e~ and 
Z —7- ji^ jjL^ are well treated with on-shell approximation. 

By using parton distribution function of gluon /g(x, s), we give the formula of the total 
cross section of pp — )• ZZ: 



_model / „\ 



f dTa^'^i''h^zz{rs)L{r,s), (71) 

JO 



where we define 



L{t,s) := dyfg{^ey,Ts)fg{V^e-y,Ts). (72) 



The concrete form of cf^^^zz been shown in Eqs. (55)-(62). The invariant mass of ZZ 
is represented as = s = ts. Then, the differential cross section of pp — t- ZZ — )• 4^ as a 

function of Mzz is written as 

''"°^!l;f '''' . ^ , X X ^ - (73, 

where Br(Z — )■ 2£) := Br(Z — t- e~^e~) = Br(Z — t- fi~^fi~) = 0.034 is the branching ratio. 



4 Numerical results 

We use the top quark mass mt = 172 GeV and the LO QCD running coupling. We apply 
the CTEQ5 LO PDF of gluon |63j . In our analysis, we consider only KK-scales large enough 
not to open up the Higgs decay channel into a pair of KK-top, Z, and W particles. When 
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Geometry 


Allowed region of KK indices 


rVZs or 


(m,n) = (1,0), (2,0), (3,0), (4,0), (5,0), 

(0,1), (1,1), (2,1), (3,1), (4,1), (5,1), 
(-1,2), (0,2), (1,2), (2,2), (3,2), (4,2), 
(-2,3), (-1,3), (0,3), (1,3), (2,3), (3,3), (4,3), 
(-3,4), (-2,4), (-1,4), (0,4), (1,4), (2,4), (3,4), 
(-1,5), (0,5), (1,5), 
(1,-1), (2,-1), (3,-1), (4,-1), (5,-1), 
(2,-2), (3,-2), (4,-2), 
^3^ —3), (4, —3), 




(m,n) = (1,0), (2,0), (3,0), (4,0), (5,0), 

(0,1), (1,1), (2,1), (3,1), (4,1), (5,1), 
(0,2), (1,2), (2,2), (3,2), (4,2), 
(0,3), (1,3), (2,3), (3,3), (4,3), 
(0,4), (1,4), (2,4), (3,4), 
(0,5), (1,5), 




(m,n) = (1,0), (2,0), (3,0), (4,0), (5,0), 

(1.1) , (2,1), (3,1), (4,1), (5,1), 

(1.2) , (2,2), (3,2), (4,2), 

(1.3) , (2,3), (3,3), (4,3), 

(1.4) , (2,4), (3,4), 

(1.5) , 


SVZ2 


j = 1 ~ 9 (n(j) = j + 1 for j: even or n(j) = j for j: odd) 


PS or 52 


j = 1 ~ 9 (2j + 1 degenerated states) 



Table 4: The region of KK summation for the maximum UV cutoff. 



the decay channel into two KK particles opens up, the Higgs resonance tends to become too 
broad and hard to be seen at Tevatron and LHC. We list in Table |4] the KK modes that 
satisfy the maximum cutoff criterion given in Table [3j 

4.1 Tevatron 

We evaluate the H — )• ZZ — )• 4^ event number per each 25GeV bin of Mzz^ expected at 
Tevatron with an integrated luminosity 8fb~^ at y/s = 1.96 TeV in the above mentioned 
UED models. We show the results for the Higgs mass Mh = 330 GeV. We consider the KK 
scales Mkk = 200, 400, and 800 GeV, except for the Dirichlet Higgs (DH) model which does 
not have a zero-mode Higgs. For the DH model, we take first-KK Higgs as the Higgs field, 
that is, Mkk = Mh = 330 GeV0 The results are shown in Fig. |4} We see that the event 
number is enhanced in all the UED models from the SM one. In particular, the 6D Projective 
Sphere model can give a large enhancement in the Higgs production by a factor as large as 



hundred compared to the SM when the KK-scale is low at 200 GeV, ^'^ Two ZZ — )• M events 

^ Note that the couphng of the Dirichlet Higgs to the SM modes is decreased by a factor 0.9 compared to 
the SM Higgs. 

Even if the one-loop H ^ gg decay rate is enhanced by the factor 100 from the SM one, it is still 
subdominant compared with the tree-level H — > WW [32] and we neglect its contribution to the total decay 
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Figure 4: The H — t- ZZ — )■ 41 event number per each 25 GeV bin of Mzz for Mh = 330 GeV 
expected at Tevatron with an integrated luminosity 8it~^ at y/s = 1.96 TeV. The grey, black, 
and magenta lines represent the expected event number with Mkk = 200, 400, and 800 GeV, 
respectively. For 6D UED models, we consider dependency on the UV cutoff, whose range is 
from minimum (lower side of band) to maximum (upper sider of band) given in Table p| 



are observed in a 300-350 GeV bin at CDF with the integrated luminosity 4.8 fb"^ [M] and 
two such events are observed in 325-375 GeV bins at DO with 6.4 fb~^ |65j . Recently two more 
events are reported around 330 GeV However, the above large cross section to explain 
the Tevatron data is found to be inconsistent with the LHC data [671 EH] . 



4.2 LHC 

We plot the event number of — t- ZZ — )■ M for the Higgs mass Mh = 500 GeV at the 
LHC with an integrated luminosity 10fb~^ at = 7 TeV in Fig. [sj When \/s = 7 TeV, 
we have checked that we cannot see sizable number of events for all the UED models with 
Mff = 700 GeV even for an integrated luminosity 10fb~^, expected by the end of 2012 
after which the LHC is planned to be (shut-down for a year and then) upgraded to 13- 

width of the Higgs. 
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Figure 5: The H —?- ZZ — > Ai event number per each 25 GeV bin of Mzz for Mh = 500 GeV, 
expected at LHC with an integrated luminosity 10 fb^^ at ^/s = 7 TeV. The black, magenta, 
and cyan lines represent the expected event number with Mkk = 400, 800, and 1200 GeV, 
respectively. Dependence on the 6D UV cutoff scale is shown the same as in Fig. [4| 



14 TeV. Therefore, we show corresponding results for Mh = 700 GeV at ^/s = 14 TeV with 
an integrated luminosity 10fb~^ in Fig. [6| We have chosen several KK scales: 400, 800, 
1200 GeV. We show our plots in logarithmic scale so that one can easily see the results for 
different luminosities by simply shifting them upward/downward. 

In Fig.jsjwith an integrated luminosity 10fb~^, we can see in total a few events in 5D UED 
models (mUED and DH) and C'(l)-C'(10) events per each 25 GeV bin for 6D UED models. 
Note that even if we can only see at best in total a few events for 5D UED by the end of 
2012, this ZZ — )• M channel is virtually background free at 500 GeV and the result would be 
still significant. 

In Fig. [6] for Mh = 700 GeV, we have plotted the results for the upgraded energy ^/s = 
14 TeV. We see that even with the integrated luminosity 10fb~^, 6D UED models can have 
a few events per each 25 GeV bin if the KK scale is relatively low Mkk = 400 GeV. The 6D 
r^/Z2, RP^, S'^ and PS models can have in total few events for higher KK mass Mkk = 
800 GeV. When the integrated luminosity adds up to 100 fb~^, we can see a few events per 
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Figure 6: The H ZZ Ai event number per each 25 GeV bin of Mzz for Mh = 700 GeV, 
expected at LHC with an integrated luminosity 10fb~^ at ^/s = 14 TeV. The black, magenta, 
and cyan lines represent the expected event number with Mkk = 400, 800, and 1200 GeV, 
respectively. Dependence on the 6D UV cutoff scale is shown the same as in Fig. [4| 



each bin for 5D mUED and DH models, and even for the SM (though the SM itself cannot 
satisfy the electroweak constraints on the S, T parameters, contrary to the UED models). We 
see that the Dirichlet Higgs model has slightly smaller cross section than the Standard Model. 
This is because the KK scale is fixed to be large, 700 GeV, and hence the enhancement of the 
KK top loop is small, while the Yukawa coupling of the (first KK) Higgs to the top quark is 
decreased by the factor 2\/2/tt ~ 0.9. 

Let us emphasize that the enhancement of Higgs production gg ^ H in UED does not 
depend on the details of the model such as the mass structure at the orbifold fixed points. 
Parameter dependence is only on the Higgs mass and the KK scale. In this sense, this 
Higgs channel signal is complementary to the direct search of the KK modes decaying into 
LKP [m [28] , which is nice because of the directness but is dependent on the details of KK 
mass splitting from the boundary terms. 
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5 Summary and discussions 



We have presented a review on the known 5D and 6D UED models focusing on the relevant 
part to the gluon fusion process. We have explained our computation of the gluon fusion 
process including the KK top-quark loops, which is new for the T^/Z2, T'^/{Z2 x Zg), RP^, 
and S'^ UED models. For 6D UED models, we have shown an NDA analysis of the highest 
possible UV cutoff scale in the 5^-based compactification, extending the analysis of Ref. [48,. 3] 
on the compactification. 

For Higgs mass Mh = 500 GeV, we can see a few (virtually background free) H — )• ZZ — )■ 
M events in 5D UED models with 10fb~^ of integrated luminosity. The 6D UED models can 
further exhibit the shape of the resonance if the KK scale is relatively low. When Higgs mass 
is as large as Mh = 700 GeV, we found no parameter region that can be seen within the 
integrated luminosity of C(10)fb~^ at ^/s = 7TeV. We have also studied the event rate for 
the upgraded energy ^/s = 14 TeV when Mh = 700 GeV. We see that the 5D and 6D UED 
models typically require 100 fb~^ and 10 fb^^ of data, respectively, in order to establish the 
existence of the resonance. As is reported in Ref. [39], the 6D UED model on the Projective 
Sphere (or S"^) shows the greatest enhancement of the Higgs production via gluon fusion 
process, among all the known UED models. 

Let us again emphasize that the presented Higgs signal of UED needs only the Higgs mass 
and KK scale as input parameters, is independent of the detailed KK mass splitting, and 
hence is unaffected by the boundary /fixed-point mass structures. Therefore, this Higgs signal 
of UED is complementary to the direct KK resonance production and to the dark matter 
signal. As the enhancement of the Higgs production via the gluon fusion process can be 
so large, recent data from the LHC |67| l68] can already significantly exclude the parameter 
space of the UED models. This analysis is presented in a separate publication [69]. We will 
also show a combined bound from the triviality and the electroweak precision constraints 
in addition to that of Ref. [69j. (The triviality bound would lower the maximum allowed 
UV cutoff scale when the Higgs mass is heavy.) In this paper, we have studied the cleanest 
possible signature H — t- ZZ — t- 4£. It is expected that a combined analysis including other 
decay channels such as H ^ WW and H — )• ZZ llvv will provide a large gain over all 
individual analyses [68j. Such a combined analysis for the UED models will be presented 
elsewhere. 
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Appendix 



A Feynman rules for Dirichlet Higgs model 



In this appendix, we describe the Feynman rules that are necessary for our computation. The 
mass terms of KK fermions are 
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mt 
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Qt 
Ut 



lin) 



(74) 



where Qt is an upper component of the quark doublet in third generation and Ut is the top 
quark singlet. Transforming each KK states by the following unitary transformation including 
chiral rotation: 
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we can obtain the ordinary diagonalized Dirac mass terms, where t^^^ and ' mass 
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top 



eigenstates of n-th KK top quarks and each mixing angle a^") is determined to be cos 2a^'"'> = 
[njli)! -^ml + m? /B?, sin2a^") = mt/ \/m1 + rfi /E?. Each KK state is twofold degenerat 
and n-th KK top mass is nit^^n) = \/rr^^^rri?jB? . The corresponding interaction terms are 
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(76) 



where (/4,5 is a dimensionlcss 4D SU{3)c coupling constant and t^EW is the 4D Higgs vacuum 
expectation value which appears after KK expansion. Q^^ ^ is massless gluon and is m-th 
KK Higgs bosons. The concrete shape of the factor of V^^n is 2-\/^/mr, whose origin is the 
non-orthonormality of mode functions in Dirichlet Higgs model. In Dirichlet Higgs model 
there is no zero mode Higgs because of choosing Dirichlet boundary condition in Higgs field. 
The first KK Higgs boson behaves like a heavy SM Higgs except that its interaction with the 
SM fields are multiplied by \/2ei. The explicit form of Feynman rules is: 
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where a is a gluon (adjoint) color index and A" are the Gell-Mann matrices. It is noted that we 
can find an interaction which is proportional to 7^ matrix at the KK Higgs Yukawa couplings 
to KK top quarks, which generates another type of contribution to the Higgs production 
through the gluon fusion process)^ This type of contribution do not exist in ordinary 5D or 
6D UED models. 



B UV cutoff based on RGE analysis in KK picture 

In this section, we show the details of the renormalization group analysis, based on the strategy 



stated in Section 3.2 that follows Ref. [48j. We note that we do not need any regularization of 
the infinite KK sum, employed in Ref. |48j , in our bottom- up approach. A higher-dimensional 
gauge theory is equivalent to the corresponding 4D theory with infinite tower of KK modes. 
It suggest that all we have to do for deriving running effect of the 4D effective gauge coupling 
at leading order is to count the d.o.f. of fields whose masses are lower compared to a reference 
energy This important information is encoded into the coefficient bj. The formula for the 
running coupling of 4D SU{N) or U{1) gauge theory is well-known as follows: 

" (81) 
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(83) 



This situation is similar to the famous fact that 7r° — > 27 decay is enhanced through chiral anomaly [TO] 
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gauge group 


SM contribution (b^^^) 


KK contribution (b^^) 


SU{3)c 


-7 


-2 


SU{2)w 


-19/6 


3/2 


U{1)y 


41/6 


27/2 



Table 5: RGE coefficients in Eq. (70). 



where is a 4D gauge coupling strength, whose group is discerned by i; C (r) is defined by 
tr [T^T^] = C (r) 5"^ for each representation of SU{N) group and the specific value is 1 /2 for 
the fundamental representation; C2 (Adjoint) is the quadratic Casimir operator for the adjoint 
representation, whose value is A^; Yf (Yh) shows U{1) charge of 4D Weyl fermion (Higgs). The 



contribution from KK vector bosons is included in the first term in Eq. (82). We note that a 
6D gauge field has two extra dimensional components and that a 4D adjoint scalar is left as 
a physical mode after the KK decomposition. (When counting the number of adjoint scalar 



degrees of freedom in the running (81), it is two rather than one in our treatment.) 



We can obtain the solution of Eq. (81 ) by the integration of the both sides over the region 
[mz,fj\- The essential point is that the coefficient bj becomes altered when the reference 
energy /i cross a threshold and the number of the effectively massless d.o.f. of fields under 
the ^ changes. Neglecting the KK mass splitting from the electroweak symmetry breaking, 
the threshold correction from l-ih KK particles arrises simultaneously when ^ exceeds the 
value of /-th KK mass. This approximation simplifies the calculation to estimate the effect 
from /-th KK particles all at once. The values of bj in the zero modes (b^) and /-th KK 
modes (bf^) from the bulk SM matter contents are summarized in Table sjp" 

This way, the shape of a4i(//) is determined as 



b^^ u b^° ^ u 
I "HO 



«4"/(/^) = «4"/("^z) - ^ In ^ - ^ ^ In (84) 

2tt mz 27r ^ mi 



where m^^) shows /-th KK mass and the upper bound of the summation is 

m^i) < /i. (85) 

In the case of in five dimensions, where the spectrum of KK masses is equally-spaced, 
the above calculation is executed with no difficulty. In 6D cases, by contrast, the KK mass 
spectrum is not equally-spaced and we use the following approximation. 

□ T^-case: In the T^-case, the form of (m,n)-th KK mass is m(^jn,n) = ^/m? + v? /R 
and eventually that of the ffist (lightest) KK mass is Mkk = 1/R- The exact form of 
the summation is as follows: 



El- 



, — for Km^+n"^ < ( ] . (86) 

, , Vm^ + n^KK ~ ~ \MkkJ ^ ' 

(m,n) 

We approximate the summation over m and n by integral over r and 9 in two-dimensional 



Note that we do not employ the GUT normahzation for the U{1)y couphng and the beta function. 
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polar coordinates: 
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(88) 



□ S'^-case: In the S'^-case, the form of (j, m)-th KK mass is ?7i(j,m) = + 1)/-^ 

and that of the first (lightest) KK mass is Mkk = V^/R- There are 2j + 1 numbers of 
degenerated states in each j-level. The exact form of the summation is as follows: 
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for 2 < j(j + 1) < 2 f — 



When we use the approximation: j{j + 1) ~ j^, jmax and jmin are determined as 



Jmax 



ttt-kk 



Jmin 



~ \/2. 



(89) 



(90) 



Let us approximate the summation over j by an integral over j: 
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Thereby we obtain the final form: 
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terms in Eq. (70), we obtain 



Combining Eq. (88) and Eq. (92), we get the final form (70). Neglecting the logarithmic 
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2vr Ml^- 



(93) 



We note that the coefficient of the quadratic term for coincide with that in Ref. [38 



obtained from a different regularization scheme. Putting Eq. (93) into the condition e(A) ~ 1, 
we get 



C{Ng + 2bDa4i{mzy 



(94) 
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Type of J 


j;eometry 


types 


T^-based 


5^-based 


SU{3)c cutoff 


no cutoff 


no cutoff 


SU{2)w cutoff 


6.9 Mkk 


8.6 Mkk 


U{1)y cutoff 


5.3 Mkk 


6.6 Mkk 


SU{3)c Landau pole 


no cutoff 


no cutoff 


SU{2)]v Landau pole 


8.9 Mkk 


11 Mkk 


U{1)y Landau pole 


5.4 Mkk 


6.8 Mkk 



Table 6: The values of the cutoff scales and the positions of the Landau poles in and S"^ 
cases. 



where we have used V2 = SttC/M^-^. Concretely, we get 



5.3 Mkk for 
6.6 Mkk for 5^ 



A<<!;_. . (95) 



: from the U{1)y cutoff. 

In addition to this analysis, we also make consideration for the Landau poles of the gauge 
interactions. If the value of the energy where a Landau pole emerges is smaller than that of 
the cutoff which we have discussed before, we should treat the position of the Landau pole: 
Q;Jj^(ALandau) = 0, which is easily obtained with leading order approximation as 

■''-Landau fiD / \ ' V'^^.' 

Chia4i{mz) 

as a cutoff scale instead. The concrete forms of each value are show in Table [H 

In the analysis above, we have taken values of Ng as 3, 2 and 1 in each case of SU{3)c, 
SU{2)\Y and [/(l)y, respectively, and have employed the values 

«c/(i)y("^^) = 97.9, 
'^SU(2)J^z) = 29 A, (97) 

at mz = 91.1 GeV. We do not consider a TeV-scale gauge coupling unification condition as 
a UV cutoff in this paper. 

In the both and cases, the most stringent bounds come from the U{1)y cutoff 
scales, which restrict the effective range of the perturbation the most severely. It is natural 
that the scale emerging the U{1)y Landau pole is near the upper limit of the perturbativity 
but a little bit higher. 

C Breit-Wignar formula 

We review how the Breit-Wigner formula emerges from the resummation of the one-particle 
irreducible (IPI) Higgs two-point function, in order to be careful of possible systematic errors 
when the width becomes broad. In this paper, we assume that the Higgs mass is larger than 
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twice the W KK mass so that it is sufficient to hmit ourselves to the SM case when discussing 
the Higgs total width. It is straightforward to extend the result for the UED when one wants 
to take the KK loops into account. 

In the SM, the Higgs production cross section via the gluon fusion process is obtained as 
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The expression (98) is obtained in the limit of the vanishing decay width Tu — t- 0. We may 

(102) 



introduce a narrow width in Eq. ( 100 ) by the Breit-Wigner type replacement 
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M^)' + A2 

When we want to reproduce the delta function, A in Eq. ( 102[ ) cannot depend on s, otherwise 
we cannot get the correct normalization: J ds6{s — Mfj) = 1. One should then perform the 
replacement of the delta function ( 102 ) as 

A ^ MhTh{Mh) (104) 

in Eqs. (98) and (100). In literature, see e.g. |32j, A in Eq. (103) is sometimes replaced as 
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Mh 
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(105) 



Instead of the truncation (100), we have already obtained the full gg ^ H ^ ZZ cross 
section (50) by the naive Breit-Wignar type replacement A — )■ itih^h in and only in the 
denominator of the Higgs propagator 

(106) 



Q^-Ml + iA' 



Hereafter, let us see that this treatment gives sufficiently good approximation to the full 
result (pel). 
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C.l Resummed propagator 

First let us review how the resummed propagator is obtained in the SM. We write the bare 
Higgs mass and field in terms of the renormalized ones and the counter terms 



The resummed bare propagator reads 

Db = 



Hb 



ZrH, 



Zh = 1 + 5Zh. 



Q^-MI + Uh{Q^ 
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(108) 



where Mb and IIh{Q^) are the bare mass and the IPI two-point function, respectively, both 
of which contain ultraviolet (UV) divergences, and := —q^ for a Higgs four-momentum q. 
The renormalized propagator then becomes 
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Zh Q^-M% + ±HiQ^y 



where 



±h{Q^) := ^HiQ^) - Zh 6Mjj + 6Zh (Q' - Mfj) 



(109) 



[110) 



is the renormalized (finite) IPI two-point function, with T,b{Q^) '■= Zb^h{Q'^) being the 
IPX two-point function that is given in terms of the renormalized fields but still contains UV 
divergences. Note that it is sufficient to consider the case > for our purpose since we 
have only s-channel Higgs propagator, though can be negative when the Higgs is virtual, 
e.g. when exchanged in t-channel. The Higgs two point function in the SM is given by [7l] 
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where the loop functions are 
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For concreteness we write down 
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I{k) = 2 + In - + ztt + 0{k In k) 



(k< 1), 

(At > 1). 



C.2 On-shell scheme renormalization 

In the on-shell scheme, we put the following renormalization conditions: 



0, 



0, 
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which gives 
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We can compare this result with the tree-level Higgs decay width 
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Figure 7: Parton level cross section for the process gg ^ H ^ ZZ with rriH = 300 GeV 
(left) and 700 GeV (right). Inclusion of full two point function (126) (blue solid line), the full 
one ([50j) with A Mh^h (red dashed), and the truncated one (l03[ ) with A sTh/Mh 
(yellow dot-dashed) are shown. In the right figure, the full one (50) with A — t- sTh/Mh (thin 
dot-dashed) and the truncated one ( 103 ) with A — )• Mh^h (thin dotted) are also drawn. 



We see from Eq. ( 123 ) that the leading term for the Q ~ Mh ^ 2mt limit is not proportional 
to for the W and Z contributions and therefore the replacement (105), namely 



Dr 



(125) 



does not give a good fit. 

When computing the full cross section for the process gg ^ H ^ ZZ ^ we may employ 



the resummed propagator (109). Neglecting the contributions from box diagrams, we get the 
cross section 



2567r2 V^EW / 
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(126) 



In Fig. we show the results for various replacement. We see that any replacement 
suffices when Higgs mass is not very large, Mh = 300 GeV (left). (This is the case for the 
Higgs mass above the top threshold mn > 2mt too.) For the large Higgs mass (right), the 



decay width becomes larger and we see that our approximation ( 50 ) with Breit-Wignar type 
replacement in the denominator A — )• Mh^h gives a good fit to the cross section with full 



two-point function (126). 



Finally, just for comparison with Eq. ( 103 ), let us present the resummed cross section ( 126 ) 
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Figure 8: Ratio of the pole mass rriH to the on-shell one Mh as a function of the on-shell 
mass (left); same for the ratio of the pole-scheme width Th to the tree- level one T^ff^ (right). 



in a rewritten form 
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C.3 Pole scheme renormalization 



Let us review the pole scheme renormalization [72]. Instead of the on-shell condition (121), 



the pole scheme renormalization condition fixes the pole of the propagator ( 108 ) at which 
we parametrize by two real constants niH and Th as: Q"^ 



m 



H 



iniH^H, namely, 



Q2-M| + nH(Q2) = o. 



(128) 



As the pole position of the bare propagator ( |108 ) is the same as that of the renormalized 
one (109), we see that the on-shell renormalized two-point function satisfies 



Q2 _ Mi + ±HiQ^) = 0, 



that is, 



Mjj = ml + RetH{Q^) 



mu^H = ImSH(<3^) 



(129) 



(130) 



We see that the real and imaginary parts determine the pole scheme mass mn and decay 
rate Th as functions of the on-shell scheme mass Mh- Note that the renormalized two-point 
function Tih has implicit dependence on the on-shell mass Mh- 

In Figjsj we plot the ratio of the pole to on-shell mass vtihIMh (left) and the pole-scheme 
width Vh to the tree- level width ( 124[ ) (right) as functions of the on-shell mass Mh^ computed 
within the SM. 
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It is known that the on-shell scheme mass is gauge dependent at the NNLO, see references 
in [72]. In contrast, the pole position of the amphtude is a gauge independent physical notion. 
Therefore, in principle we should utilize the pole scheme mass and width. However, the on- 
shell mass and the pole mass are identical at consistent 1-loop order (which is the highest 
order considered here), and they start differing only at 2-loop order and above. We see 
from Fig. [s] that the both schemes agree within 1% accuracy in our approximated treatment 
neglecting non-resonant box contributions, which are the same order as 0(a)-corrections to 



the resonant ones,^'^ We safely utilize the on-shell mass and the tree-level width (124) even 
though the Higgs decay width becomes as large as 180 GeV when Mh = 700 GeV. 

Note that our treatment to include the decay width in the numerical calculation (|50|) 



with (54) corresponds to the (gauge invariant) complex-mass scheme, see e.g. Refs. |73ji74j, 
where the real Higgs mass is replaced by a complex value everywhere in the amplitude, which 
in our case leads to a replacement of the Higgs mass in its propagator only. 
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